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1 Introduction 



Here we present formulae for derivatives of associated Legendre functions (hereafter referred to as Legendre functions) 
of the first kind P„{z) and the second kind Q£(z), with respect to their parameters, namely the degree v and the 
order fi. Some formulae relating to these derivatives have been previously noted [TO] and also there has been recent 
work in this area [2l [T5j [13l [14l [12] with Brychkov (2009) [3] giving a recent reference covering the regime for argument 
z G [—1,1]. In this paper we cover parameter derivatives of Legendre functions for argument \z\ > 1. 

The strategy applied in this paper is to incorporate derivatives with respect to order evaluated at integer orders for 
modified Bessel functions (see [JJ|TJ[T0]) to compute derivatives with respect to the degree and the order of Legendre 
functions. Below we apply these results through certain integral representations of Legendre functions in terms of 
modified Bessel functions. 



2 A useful transformation on the complex plane 

There is a transformation over an open subset of the complex plane which is particularly useful in studying Legendre 
functions Oil]. This transformation, which is valid on a certain domain of the complex numbers, accomplishes the 
following 

cosh z <-> coth w 
coth z «-> cosh w 
sinhz «-> 1/ sinhiu, 

where cosh, coth, and sinh are the complex hyperbolic cosine, cotangent, and sine functions respectively. This trans- 
formation is accomplished using the following map 

w(z) = log coth-, (1) 

(where log is the complex natural logarithm) which is verified to be an involution over an open subset of the complex 
plane given by — tt < Im(z) < w, where one removes a branch given by Re(z) < 0, Im(z) = 0. This mapping is 7r 
periodic in the imaginary direction and is locally injective over the entire complex plane when the following set of 
complex numbers are removed {z : z — inn, n 6 Z} U {z : Im(z) = 2nn and Re(z) < 0, n G Z, z € C}. 

This transformation is particularly useful for certain Legendre functions which have natural domain given by the 
real interval (1, oo), such as toroidal harmonics [HllS] (and for other Legendre functions which one might encounter in 
potential theory), Legendre functions of the first and second kind with odd-half-integer degree and integer order. The 
real argument of these Legendre functions naturally occur in [l,oo), and these are the simultaneous ranges of both 
the real hyperbolic cosine and cotangent functions. 

One application of this map occurs with the Whipple transformation of Legendre functions [161 [7] under index 
(degree and order) interchange. See for instance, eqs. (8.2.7) and (8.2.8) in [JJ, namely 

-"-^KVF^TJ ~ Vtt r^ + M + i) QA h [ ) 

which are valid for Re z > and for all complex v and fi, except where the functions are not defined. T(z) is the 
Gamma function (see PQ). 

3 Parameter derivative formulas from K v {t) 

Starting with Gradshteyn & Ryzhik (2007) eq. (6.628.7) [8] we have 

£ e-**K v {t)t a - 1/2 dt = ^lv{a-v+ 1 ^T{a + u+ l ^{z"- lp /2 P^ l/2 {z) 

- r (a - , + i) - 1) e-Q^ 1/2 , (3) 
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where K v (t) is a modified Bessel function of the second kind with order v, and the two equalities are established 
through the Whipple transformation, eq. @. 

We would like to generate an analytical expression for the derivative of the Legendre function of the second kind 
with respect to its order, evaluated at integer orders. In order to do this our strategy is to solve the above integral 
expression for the Legendre function of the second kind, differentiate with respect to the order, evaluate at integer 
orders, and take advantage of the corresponding formula for differentiation with respect to order for modified Bessel 
functions of the second kind (see [31 (H [H [10] ) . Using the expression for the Legendre function of the second kind in 
eq. (|3|). we solve for Q t ^_ 1 / 2 ( z ) an d re-express using the map in eq. ([1]). This gives us the following expression 



(*» - 1) 



-I//2-1/4 



,lTTfJ, pOO 



exp 



Differentiating with respect to the order fj, and evaluating at ji = ±m, where m £ Nq = N U {0} yields 



fi—^zm 



= ^ i) 



-J//2— 1/4 







(*» - 1) 



-v/2-l/A 



(-1) 



m poo 



r(z/ T m+i) 
The derivative from the first term is given as 
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r \V =p TO 
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where the Digamma function ip(z) is defined in terms of the derivative of the Gamma function with respect to its 
argument z through 

T'(z) = T(z)i>(z). 
The derivative in the second integral (see [31 El [IJ [10]) is given by 

m— 1 



d_ 



±to! 



1 



j-k — m 



k\{m - k) 2 k - m + 



~K k {t) 



(4) 



(see for instance eq. (1.14.2.2) in Brychkov (2008) [3]). Substituting these expressions for the derivatives into the two 
integrals and using the map in eq. fl} to re-evaluate these integrals in terms of Legendre functions gives the following 
general expression for the derivative of the Legendre function of the second kind with respect to its order evaluated 
at integer orders as 



T[y-m+ f) 



fj,— ±m 
k-m ( -2 i\( fe - m )/ 2 



I a f ' ( U^f TO + — 



Q™i /2 (z) 



' ^ k\(m- k)2 k - m+1 



~Qv+k-m-l/2( Z )- 



For fj, — there is no contribution from the sum and the result is 
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in + tp \ v + 



Qv-l/2{z), 



which agrees with that given in Magnus, Oberhettinger, and Soni (1966) [10] in §4.4.3. We are now able to obtain 
formulas for non-zero values of ji such as for fi = — 1 



v 
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in + tp ( v + - 



Ql- 1/2 (z) + {z 2 - 1) 1/2 Q v - m (z) 
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or for /i = +1 



ITT + 1p \ V — — 



or for other integer values of \x. 

If we start with the expression for the Legendre function of the first kind in eq. ([3]) and solve for (2) we nave 



(z 2 - 1) 



M/2 



e- zt K u (t)t^ 1/2 dt. 



- 1 ' 2y "' V 7r r (// — z/ + 1) r (// + + 1) j 

Differentiating with respect to degree v and evaluating at v — ±n, where ngNo yields 



(5) 



d 



V— ±71 



M/2 



i 



^ r (// - + 1) r (fj, + z/ + \ 



/"OC 

Jo 



dv 



K u {t) 



dt. 



v—^sm 



The derivative from the first term is given as 
d 1 



dv T (/x - i/ + |) r (/i + i/ + |) 



_ V (m =F " + I) - ip (/x ± n + |) 



Substituting this expression for the derivative and that given in eq. ((4J yields the following general expression for the 
derivative of the Legendre function of the first kind with respect to its degree evaluated at odd-half-integer degrees as 



± 



r(/>- 



v— ±n 
n-1 

2/ k=0 



r (/x - n + 2A; + |) (z 2 - l) 



(n-fc)/2 



fe!(n- /c)2 fc - 



n+l 
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(z). 



If one makes a global replacement, —fx t— ► /x, using the properties of Gamma and Digamma functions, this result 
reduces to 



± 



c) 



v— ±n 
n-1 
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Note that by using the recurrence relation for Digamma functions (see §1.2 in [TO] 

#z + l) = + 



?/> ( ix + n + i J - V ( A* - « + \ ) = 2^^] 



we can establish 
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For v — there is no contribution from the sum and the result is 
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which agrees with that given in Magnus, Oberhettinger, and Soni (1966) [TUj in §4.4.3. We are now able to obtain 
formulas for non-zero values of v such as for v = ±1 
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i/=±i 



2/iPf /2 (z) + (z 2 



i) 1/2 ^r; 2 w, 



or for other integer values of v. 

Note that this method does not seem amenable to computing derivatives with respect to the degree of Legendre 
functions of the form P^(z) evaluated at integer degrees, since shifting the degree by +1/2 in eq. (JSJ) converts the 
modified Bessel function of the second kind to a form like K^ + i/ 2 (t) and the derivative with respect to order of this 
Bessel function (see [H [TJ [10]) is not of a form which is easily integrated. 

4 Parameter derivative formulas from I u (t) 

Starting with Gradshteyn & Ryzhik (2007) eq. (6.624.5) [3], see also Prudnikov et. al. (1988) [UJ, we have 



%(t)t Q - 1/2 di 
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where I v {t) is a modified Bessel function of the first kind with order v 1 and the two equalities are established through 
the Whipple transformation, eq. @. We will use this particular integral representation of Legendre functions to 
compute certain derivatives of the Legendre functions with respect to the degree and order. 

We start with the integral representation of the Legendre function of the second kind in eq. ©. Differentiating 
with respect to the degree v and evaluating at v = ±n, where n £ No, one obtains 



u—zen 
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The derivative of the modified Bessel function of the first kind in eq. ([7]) (see [21 IU [TJ [10]) is given by 
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(see for instance eq. (1.13.2.1) in Brychkov (2008) [3]). Inserting eq. ([5]) into eq. Q and using eqs. ^ and ([5]), we 
obtain the following general expression for the derivative of the Legendre function of the second kind with respect to 
its degree evaluated at odd-half-integer degrees as 
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For v = there is no contribution from the sum and the result is 
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which agrees with that given in Magnus, Oberhettinger, and Soni (1966) [10] in §4.4.3. We are now able to obtain 
formulas for non-zero values of v such as for v = ±1 



v=±\ 



= -d~^T L-i) (z 2 - 1)- 1/4 Ql_ 1/2 



±(z 2 -l) 1/2 Q^Uz), 
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or for other integer values of v. 

We can see that this method does not seem amenable to computing derivatives with respect to the degree of 
Legendre functions of the form Q$(z) evaluated at integer degrees, since shifting the degree by +1/2 in eq. ([7]) 
converts the modified Bessel function of the first kind to a form like I v +i/ 2 (t) and the derivative with respect to order 
of this Bessel function (see [H [TJ [10]) is not of a form which is easily integrated. 

Finally, we obtain a formula for the derivative with respect to the order for the Legendre function of the first kind 
evaluated at integer orders. In order to do this we use the integral expression for the Legendre function of the first 
kind given in eq. (|6|) and the map given in eq. ([1]) to convert to the appropriate argument. Now use the negative order 
condition for Legendre functions of the first kind (sec for example eq. (22) in iTj) to convert to a positive order. 

Differentiating both sides of the resulting expression with respect to the order fi and evaluating at \x = ±m, where 
m G N yields 
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The derivative of the reciprocal of the Gamma function reduces to ip (v =p m + 1/2) /F (y =F m + 1/2). The derivative 
with respect to order for the modified Bessel function of the first kind is given in eq. ([5]). The integrals are easily 
obtained by applying the map given by eq. |T]) as necessary to eqs. ^ and Hence by also using standard properties 
of Legendre, Gamma, and Digamma functions we obtain the following compact form 
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For fi = there is no contribution from the sum and the result is 
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which agrees with that given in Magnus, Oberhettinger, and Soni (1966) [10] in §4.4.3. We are now able to obtain 
formulas for non-zero values of ji such as for fi = — 1 



dfl P u-l/2( Z ) 



= Ql-i /2 (z) + V (v + \ ) P^ 1/2 (z) + (z 2 iy 1/2 P v - m {z) 
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or for other integer values of /i. 



5 Acknowledgements 

I owe much thanks to the following people: Prof. Radoslaw Szmytkowski suggested that I submit this work for 
publication and participated in several important discussions. I had many valuable discussions with Prof. Ernie 



6 



Kalnins, who supported me when I was writing this article at the University of Waikato. Dr Yury Brychkov engaged 
me in valuable discussions. At the University of Auckland, I had some valuable conversations with Dr Garry J. Tee 
who carefully proofread this article, and Dr Tom ter Elst discussed matters relating to the map in §2. And I thank 
the referee for helpful remarks. 

References 

[1] M. Abramowitz and I. A. Stegun. Handbook of mathematical functions with formulas, graphs, and mathematical 
tables, volume 55 of National Bureau of Standards Applied Mathematics Series. For sale by the Superintendent 
of Documents, U.S. Government Printing Office, Washington, D.C., 1964. 

[2] Yu. A. Brychkov. Handbook of Special Functions: Derivatives, Integrals, Series and Other Formulas. Chapman 
& Hall/CRC Press, Boca Raton-London-New York, 2008. 

[3] Yu. A. Brychkov. On the derivatives of the Legendre functions Pjf{z) and Q£(z) with respect to /i and v. Integral 
Transforms Spec. Fund., 21:DOI-10. 1080/10652460903069660, 2010. 

[4] Yu. A. Brychkov and K. O. Geddes. On the derivatives of the Bessel and Struve functions with respect to the 
order. Integral Transforms Spec. Fund., 16(3):187-198, 2005. 

[5] H. S. Cohl, A. R. P. Rau, J. E. Tohlinc, D. A. Browne, J. E. Cazes, and E. I. Barnes. Useful alternative to 
the multipole expansion of 1/r potentials. Physical Review A: Atomic and Molecular Physics and Dynamics, 
64(5):052509, Oct 2001. 

[6] H. S. Cohl and J. E. Tohlinc. A Compact Cylindrical Green's Function Expansion for the Solution of Potential 
Problems. The Astrophysical Journal, 527:86-101, December 1999. 

[7] H. S. Cohl, J. E. Tohlinc, A. R. P. Rau, and H. M. Srivastava. Developments in determining the gravitational 
potential using toroidal functions. Astronomische Nachrichten, 321(5/6):363-372, 2000. 

[8] I. S. Gradshteyn and I. M. Ryzhik. Table of integrals, series, and products. Elsevier/ Academic Press, Amsterdam, 
seventh edition, 2007. Translated from the Russian, Translation edited and with a preface by Alan Jeffrey and 
Daniel Zwillinger, With one CD-ROM (Windows, Macintosh and UNIX). 

[9] E. W. Hobson. The theory of spherical and ellipsoidal harmonics. Chelsea Publishing Company, New York, 1955. 

[10] W. Magnus, F. Obcrhcttinger, and R. P. Soni. Formulas and theorems for the special functions of mathematical 
physics. Third enlarged edition. Die Grundlchren der mathcmatischen Wissenschaften, Band 52. Springer- Verlag 
New York, Inc., New York, 1966. 

[11] A. P. Prudnikov, Yu. A. Brychkov, and O. I. Marichev. Integrals and series. Vol. 2. Gordon & Breach Science 
Publishers, New York, second edition, 1988. Special functions, Translated from the Russian by N. M. Queen. 

[12] R. Szmytkowski. Addendum to: "On the derivative of the Legendre function of the first kind with respect to its 
degree". J. Phys. A, 40(49): 14887-14891, 2007. 

[13] R. Szmytkowski. A note on parameter derivatives of classical orthogonal polynomials. ArXiv e-prints, January 
2009. 

[14] R. Szmytkowski. On the derivative of the associated Legendre function of the first kind of integer degree with 
respect to its order (with applications to the construction of the associated Legendre function of the second kind 
of integer degree and order). Journal of Mathematical Chemistry, 46(l):231-260, June 2009. 

[15] R. Szmytkowski. On the derivative of the associated Legendre function of the first kind of integer order with 
respect to its degree. ArXiv e-prints, July 2009. 

[16] F. J. W. Whipple. A symmetrical relation between Legendre's functions with parameters cosh a and cotha. 
Proceedings of the London Mathematical Society, 16:301-314, 1917. 

Note added in proof. — It has recently come to our attention that the argument domain of applicability of the 
formulas for the derivatives of the Legendre functions that we have presented in this paper are actually valid in the 
complex domain z e C \ [—1, 1] D {z : \z\ > 1, z G C}. 



7 



